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Abstract 

Superscaling analyses of few-GeV inclusive electron scattering from nuclei are extended to include 
not only quasielastic processes, but now also into the region where A-excitation dominates. It is 
shown that, with reasonable assumptions about the basic nuclear scaling function extracted from 
data and information from other studies of the relative roles played by correlation and MEG 
effects, the residual strength in the resonance region can be accounted for through an extended 
scaling analysis. One observes scaling upon assuming that the elementary cross section by which 
one divides the residual to obtain a new scaling function is dominated by the ^ A transition 
and employing a new scaling variable which is suited to the resonance region. This yields a 
good representation of the electromagnetic response in both the quasielastic and A regions. The 
scaling approach is then inverted and predictions are made for charge-changing neutrino reactions 
at energies of a few GeV, with focus placed on nuclei which are relevant for neutrino oscillation 
measurements. For this a relativistic treatment of the required weak interaction vector and axial- 
vector currents for both quasielastic and A-excitation processes is presented. 

PACS numbers: 25.30.Pt, 23.40.Bw, 24.10.Jv 

Keywords: Charge-changing neutrino reaetions. scaling, quasielastic and Delta excitation 
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I. INTRODUCTION 



In recent studies of inclusive electron scattering at intermediate to high energies from 
nuclei we have explored various aspects of scaling and superscaling . The 

general procedure used in such analyses is to divide the experimental inclusive cross section 
by an appropriate single-nucleon cross section, having contributions from Z protons and 
neutrons with their corresponding electromagnetic form factors, to obtain a reduced cross 
section. Here the word "appropriate" entails two things. First, the usual analysis in the 
region of the quasielastic (QE) peak assumes that the dominant process is elastic scattering 
from nucleons in the nuclear ground state followed by quasi-free ejection of the nucleons 
from the nucleus, and hence the appropriate single-nucleon form factors are the elastic 
ones. Second, the nucleons in the nuclear ground-state are moving ("Fermi motion") and 
accordingly the single-nucleon cross section used must take this into account. Once one has 
the reduced cross section it can be plotted versus one or more appropriately chosen variables: 
if the results do not depend on some of these variables and a universal behavior is found 
one says that the results scale. Specifically, when the reduced cross section is plotted versus 
a well-chosen scaling variable (see below) and no dependence on the momentum transfer q 
is observed one says that one has scaling of the first kind. When no dependence occurs on 
the momentum scale that characterizes specific nuclei (essentially the Fermi momentum kp 
of a given nuclear species) one says that one has scaling of the second kind. If both types 
of scaling behavior are found one says that superscaling occurs. 

At sufficiently high energies we have seen both types of scaling behavior. For specific 
nuclei one observes quite good first-kind scaling at excitation energies below the QE peak, 
namely, in the so-called scaling region. This is the familiar y-scaling behavior. At energies 
above the peak, where nucleon resonances (especially the A) are important, this type of 
scaling is broken for the total reduced cross section. On the other hand, from what data 
we have where longitudinal-transverse separations have been made, we know that these 
scaling violations apparently reside in the transverse response, but not in the longitudinal. 
The latter appears to superscale. In fact, this is not unexpected, since we know that there 
are contributions which do not scale arising from meson-exchange currents (MEG) plus the 



3 



correlation effects^ required by gauge invariance which must be considered together with the 
MEC BflQ and from inelastic scattering from the nucleons p]. As discussed below, it is 
important to observe that MEC and inelastic contributions are predominantly transverse in 
the kinematic regions of interest in the present work. Note that scaling of the second kind 
works very well in the scaling region and, even in the resonance region, is only violated at 
roughly the 20% level. 

Using these recent studies as a basis we now extend our analysis to encompass both the 
QE and A regions. Our approach is the following: taking all of the high-quality data for two 
specific nuclei of relevance for our later discussions, namely for ^^C and ^^O, we proceed as 
outlined above and begin by performing a comprehensive scaling analysis in the QE region. 
Using our knowledge of the experimentally determined longitudinal superscaling function 
as the starting point, we work backwards and predict the transverse response one would 
obtain strictly from the contributions that are present in that function. In other words, we 
reconstruct the part of the transverse cross section that does not have either MEC effects 
or inelasticity built in. The net inclusive cross section so obtained is then, in effect, the 
QE contribution, except for corrections arising from MEC and their associated correlations. 
The next step is to subtract this from the data. What is left should now be dominated by 
the inelasticity, and, in particular, when not too far above the QE peak region, we expect 
the A resonance to be the most important contribution. 

From other studies, we expect that a similar procedure can now be followed for the 
subtracted results. We again reduce the left-over cross section by dividing by the appropriate 
single-nucleon cross section, now for the ^ A transition, and display the result versus a 
new scaling variable in which the kinematics of resonance electro-production are respected. 
As discussed later, the results scale quite well, suggesting that this procedure has indeed 
identified the dominant contributions not only in the QE region, but also in the A region. 
We check our analysis by assembling all of the pieces obtained via the scaling procedures 
to produce a total inclusive cross section which can be compared with data. Overall the 
answers are very encouraging, and only for specific kinematics do we see deviations as large 

^ Note that such correlation effects are not all that enter: in particular, even in a factorized approach of 
the type presented here there are both mean-field and short-range correlations in the initial state which 
are embodied in a nuclear spectral function and which lead to the scaling function deduced from data, as 
discussed below. 
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as 10-20%. We believe that effects from MEC and their associated correlations (which are 
not incorporated using this approach) are probably responsible for these residuals. 

Having met with success in extending the scaling and superscaling analyses from the 
scaling region, through the QE peak region and now into the A region we are in a position 
to take a step in a new direction. Since we have the scaling functions and can be sure that, 
upon being multiplied by the electromagnetic N ^ N and N ^ A single-nucleon cross 
sections, the total nuclear electromagnetic cross section is quite well reproduced, we can 
just as well multiply by the corresponding charge-changing weak interaction N ^ N and 
— > A single-nucleon cross sections to obtain predictions for neutrino reactions in nuclei 
at similar high-energy kinematics. 

Thus, the second motivation for the present investigation has been to work backwards to 
predictions for these cross sections with the goal of providing high-quality results for use in 
on-going experimental studies of neutrino oscillations at GeV energies 3, 12, 3, JJ, Q|. 
These studies are presently being pursued in the MiniBooNE and K2K/T2K experi- 
ments 111. Both of these and also the forthcoming MINOS, NOvA and MlNERvA 
mm .nit.a.ve. involve ne..n„o ene.gie. of .ve.a, GeV „he. a Ml. .,a...ic 
treatment of the neutrino-nucleus scattering is mandatory, but hard to achieve. Targets of 
hydrocarbon or water are involved in the cases of MiniBooNE and K2K/T2K, and hence 
^^C and ^^O are taken as the focus in the present work. For the others, iron and lead will 
also be considered, and in this regard we note that, to the degree that scaling of the second 
kind is reasonably well satisfied, one can focus on the nuclei where the most reliable electron 
scattering data are available and then subsequently obtain predictions for neutrino reactions 
not only for those nuclei, but for a wide range of targets. 

Any reliable calculation for neutrino scattering should first be tested against electron 
scattering data. Here, instead of using a specific model to describe inclusive electron scat- 
tering at relatively high energies in the QE and A regions, as stated above, we follow a 
different approach. Based on the scaling behavior of the electron-nucleus cross section in 
both the QE and A peaks we extract the scaling functions directly from experiment and use 
them to predict the neutrino-nucleus cross section. 

.o.».ed by .he fac. tha, wM,e .odeU,. Q H Q 
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25l ] of the nuclear dynamics in studies of high-energy inclusive lepton 



scattering is expected to be capable of getting the basic size and shape of the cross section. 
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so far it has not been capable of accounting for important details of the response. Specifically, 
such modeling has been able to provide a reasonable representation of the eA inclusive cross 
section. Typically at high energies the cross sections obtained using wide classes of models, 
including those with relativistic mean-field dynamics and RPA-type correlations included, 
are seen to be very similar to the results found using the Relativistic Fermi Gas (RFG) model 
and accordingly we will also make comparisons of results obtained using the scaling approach 
with those obtained using the RFG model. At the peaks of the QE or A responses one 
finds that the two approaches differ by about 25% (with mean-field effects this discrepancy 
is reduced to perhaps 20%); however, as we shall see later, the phenomenological scaling 
approach requires a long tail which is largely absent in most modeling. A possible reason 
for this disagreement is the absence of classes of short-range correlation effects in most of 
the relativistic modeling. It goes without saying that non-relativistic modeling is completely 
inadequate at the energies of interest in the present work. At intermediate energies (below 
those considered in this study) it is, of course, important to include effects from final-state 
interactions and RPA correlations (see 
can be significant. 

Since most semi-leptonic scattering processes at similar kinematics have much the same 
character, one should expect that failure at this level will also imply a similar level of 
disagreement for predictions made of neutrino reaction cross sections using the same types 
of modeling. Clearly, on the one hand, if all one wants is a rough estimate of neutrino 
reaction cross sections at similar energies, then existing relativistic modeling is probably 
adequate. On the other hand, if uncertainties of less than 25% are required (as, for example, 
when one wishes to see distortions in the energy distributions of the detected muons in 
[u^, fi~) reactions with nuclei caused by neutrino oscillations), then one must use existing 
models with great caution. As we shall see below, it appears that the superscaling approach 
being followed in the present work is capable of reducing the uncertainty to perhaps the 10% 
level, at least when one limits the focus to the QE region and the region up to inelasticities 
where the A contribution reaches its maximum. 

The paper is organized the following way. In Sec. |n] we begin with a brief discussion 
of kinematics, since we will be interested not only in electron neutrino induced reactions, 
where the lepton masses can safely be ignored, but also in muon neutrino induced reactions 
where the energies, although relatively high, are not high enough to safely ignore the muon 
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mass. In Sec. IIIII we present a summary of the formalism needed in studies of scaling and 
superscaling, both for the QE and A regions. There we give the results of our analysis 
of inclusive high-energy inelastic electron scattering data for carbon and oxygen, using the 
procedures outlined above, and thereby validate the scaling functions we use in the rest of 
the paper. In Sec. llVl we turn to the second theme of the paper and discuss charge-changing 
neutrino reactions with nuclei. We begin by presenting the basic formalism required in 
treating electroweak processes, followed by development of the single-nucleon responses in 
the QE region (Sec. IIV A|) and in the A region (Sec. IIVB|) . now of course with both vector 
and axial- vector N ^ N and — > A currents. Section FlV (^1 then contains a discussion of 
the formalism involved in obtaining the cross sections and response functions. Once these 
developments are in hand, we proceed in Sec. |3 with a presentation of our predictions for 
charge- changing neutrino reactions with nuclei. Finally, in Sec. IVII we summarize our work 
and present our conclusions. 

II. LEPTON SCATTERING KINEMATICS 

In this section we begin with a brief discussion of the kinematics involved in studies 
of lepton scattering from nuclei, including electron scattering and the subject of Sec. IIVI 
charge- changing neutrino reactions. We start with a general scattering problem in which 
an incident beam of leptons with 4-momentum K'^ = (e, k) scatters and a lepton with 
4-momentum K''^ = (e', k') emerges. In general one has 

e = V m? + k? (1) 
e' = + (2) 

where m and m' are the masses of the incident and outgoing leptons, respectively. Clearly, 
for electron scattering m = m' = rrie (usually, but not always, this can be taken to be zero) 
and for electron neutrino induced charge-changing neutrino reactions m = m^e = 0, whereas 
m' = rrif, (again, essentially zero). The difficult case is for muon neutrino induced charge- 
changing neutrino reactions where m = m^^ = 0, whereas m' = the last is clearly not 
negligible for the kinematics of interest in the present work. 
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As usual one has a 4- momentum transfer = {uj,q) with 



u! = e — e 
q = k-k', 



(3) 

(4) 



the energy transfer and 3-momentum transfer, respectively. The momentum transfer is 
spacelike: —Q^ = — u;^ > 0. For convenience we define an average leptonic mass 



M 



- (m2 + m'2) > 
2 V J - 



(5) 



and, given an excitation from target rest mass Mj to some final rest mass Mf > Mi (that 
is, the final hadronic rest frame total energy is W = Mf), define a sort of excitation energy 

1 



2M 



(Mj - Mf) > 0. 



(6) 



Then from energy-momentum conservation one has 

m 



2M; 



(7) 



Solving Eqs. Q and ((Tj) together one can get expressions for the scattered lepton's energy 
and 3-momentum. Defining 



^2 



Mf + 2Mie + w? + k'^ sin^ 9 

^/M,,{e-uJo) + M\ 



(8) 
(9) 



where 9 is lepton scattering angle (the angle between k and k'), it can be shown that 
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elikcosO) + {Mi + e)\ el - m'hj 



4 



(10) 
(11) 



el {Mi + e) + {k cos ^) y e| - m'^^ej 
where for the results to be real for all scattering angles the beam energy must be greater 

m'uo + (m'^ — w?)/2 



than emm, where 



Cmin =m' +lxlo + 



(12) 



Mi - m' 

Hence, for a given excitation energy ujq and for given beam energy and scattering angle, the 
quantities ei,2 can be computed and through them the final lepton's energy and 3-momentum 
are fixed. Clearly the 4-momentum transfer is then given as well. 
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III. SCALING AND SUPERSCALING 



A. Scaling in the quasielastic peak region 

In this subsection we briefly review the structure of the nuclear responses in the region of 



lat has 

urn 



the quasielastic peak (QEP). We begin with the basic Relativistic Fermi Gas model t 
been used to motivate scaling and superscaling behavior in this region of kinematics 
0]. Here a single parameter characterizes the dynamics, namely, the Fermi momentum kp- 
In the present work our goal is to use the electron scattering cross sections as input, perform 
a scaling analysis, and arrive at predictions for the charge-changing neutrino cross sections. 
Accordingly the focus is placed on kinematic regimes where the cross sections (induced by 
electrons or neutrinos) are substantial, and from past work it is known that under such 
circumstances it is a good approximation to work only to leading order in an expansion in 
rjp = kp/rajsi. Also = y^l + rjp — 1 = rip/2 is small. 

The leading-order QE responses (denoted bv subscript zero) may be written, in the non- 
Pauli-blocked domain, in the following form jj: 



i??^(K, A)o = Ao- [(1 + r) W2ir) - W^i(r)] x /rpgW (13) 

T 

/2?^(k, A)o = Ao [2Wi{r)] x fupaW, (14) 

with 

^0 = 3" = TTl"' (15) 

rriNKrip Inkp 

and Wi, W2 the structure functions for elastic scattering. As usual the proton [M = Z) and 
neutron (A/" = A^) contributions should be separately computed with the appropriate form 
factors and added together. The latter are linked to the Sachs form factors through the 
well-known relations 

{l + r)W2{r)-W^{r) = GI{t) (16) 
2W,{t) = 2tGI{t) . (17) 

As usual, we employ dimensionless variables A = a;/2mjv, k = q/2mi^ and r = |(5^|/4m^ = 
fi;^ — A^. The RFG has the following universal form for the superscaling function: 

fRPG{i') = \{i-i^')e{i-i^'), (18) 



that is, when plotted against the scahng variable ip, 

1 A-r 



(19) 



^ [l + X)r + K^^^iTr) 

a universal behavior is obtained with no dependence left either on momentum transfer 
(scaling of the first kind) or on nuclear species via kp (scaling of the second kind). 

In studies of electron scattering scaling one usually includes a small energy shift by 



replacing uj hj uj — E^hif ihi order 
tp' = (see, for example, pl y, IM 



to force the maximum of the QE response to occur for 
I). This is equivalent to taking A ^ A' = A — Xshift with 
Khift = Eshiftf^^N and correspondingly t ^ r' = — A'^ in Eq. (fTn|) . Often we shall use 
the pair of variables (k, ip') in place of the original pair (g, u) to characterize the inclusive 
scattering responses — clearly they are functionally related by the above equations. 

Using the guidance provided by the RFG, the procedure to be adopted in order to get 
the experimental scaling function F'^^{n^ip') in the QE domain is then clear: one simply 
divides the experimental QE cross section by 



(20) 



where ctm is the Mott cross section, vl, vt are the kinematic factors defined below in 
Eqs. fl74l75p and where the functions are 



qQE 

(Jrp 



2r 

T 



[ZGX + iVGlJ + 0{ril) 



-[ZGl, + NGl,n]+0{vl). 

K 



(21) 

(22) 



The factors involving k and r in 



Cqs. ()21|) and (j^^ arise partly from the Jacobian of 
the transformation from A to ^ j7| and partly from the explicit calculation leading to 
Eqs. p3|) and (fT^ . Finally, as in past discussions of scaling of the second kind, one multiplies 
F^^{k,iP') by kp to obtain the superscaling function f'^^{K,ip'). 
The nuclear response functions all have the general structure 



[R] 



QE 



(23) 



where A/" is the appropriate nucleon number. In particular, as stated above, one copy of 
this expression with proton form factors and M = Z should be added to another with 
neutron form factors and M = N for electron scattering. Here [i?]*'"" is the corresponding 
single-nucleon response. 
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In previous work P, [s], ^ we have shown that scahng works quite well for all nuclei and 
for energy loss uj < ujqep where ojqep corresponds to the maximum of the quasielastic peak 
— the scaling function F^^{k,iIj') is indeed largely independent of the momentum transfer 
as long as q is of the order 2k p or larger. Deviations from scaling, which mainly occur at 
larger energy loss, are related to contributions beyond quasielastic scattering such as those 
from meson exchange currents and A-excitation (see below). 

The scaling behavior becomes particularly clear if one studies the experimental response 
separated into its longitudinal (charge) and transverse (magnetic) pieces. The non-scaling 
contributions mentioned above mainly occur in the transverse response. Accordingly, one 
finds that the experimental longitudinal responses scale much better, and to much larger 
energy loss. The approach taken in ^, [4 1 therefore has been to use the experimental longi- 
tudinal responses to define the scaling function f^^{ip'). 

The total inclusive electron scattering response is then assumed to be composed of several 
contributions: (1) there is the entire longitudinal contribution which appears to superscale 
and to be represented by the empirical scaling function f^^{ip')'i (2) part of the transverse 
response arises from quasielastic knockout of nucleons from the nucleus and is also driven 
by the scaling function f^^{ip')] however, (3) the transverse response has additional con- 
tributions, at least from MEC effects with their associated correlations and from inelastic 
single- nucleon processes including the excitation of the A. From our past work we know that 
typically the effects under item (3) break the scaling. The contributions from MEC together 



Qy,yy,y,y,y,y,Q 



with their attendant correlations enter roughly at the 10% level 
and, as we argue later, may be less important for charge- changing neutrino reactions than 
they are for electron scattering — accordingly we shall ignore these non-scaling effects in 
the present work. In future work we hope to address this issue more directly with continued 
relativistic modeling of these contributions. Other effects can enter into the dynamics and 
invalidate the picture here (for example, at low q RPA correlation effects can modify the 
longitudinal and transverse responses in different ways because of the very different isospin 
character of these two channels); however, again, at the kinematics of interest in the present 
work these effects are thought to be relatively small. 

The largest non-scaling contribution to the transverse response is then believed to be the 
one arising from inelastic, but impulsive processes, especially via the excitation of the A for 
the kinematics of interest in the current study, and this provides the focus for the following 
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subsection. 

In thfi inter-comparison of the scaling functions for various nuclei has been performed 
in terms of the functions f^^{n^ ip') extracted from F^^{n^ ip') as discussed above. Excellent 
scaling of the second kind, i.e., scaling functions f^^{K.,il)') that closely match for different 
nuclei, was observed and, indeed, such second-kind scaling is actually significantly better 
realized than is scaling of the first kind. 

The combination of the scaling of the first and second kind — superscaling — allows 
one to determine from the data a universal scaling function /^^{ip'). The scaling function 
(and quasielastic cross section) for individual nuclei can then be recalculated once the Fermi 
momentum of the nucleus is known. 

0.8 I — ^ — ^ — ^ — ^ — I — ^ — ^ — ^ — ^ — I — ^ — ^ — ^ — ^ — I — ^ — ^ — ^ — ^ — I 




FIG. 1: Averaged experimental /'^^(V'q_e) versus iP'qe in the quasielastic 
region together with a phenomenological parameterization of the data. The 
integral of the curve has been normalized to unity. 



Reliable separations of data into their longitudinal and transverse contributions for A > 4 
are available only for a few nuclei js^; all of these response functions have been used to 
extract the "universal" quasielastic response function and to obtain a parameterization 
by a simple function. Figure ^ shows f^^i'ipQE) averaged over the nuclei employed, together 
with the corresponding fit^. Note that f'^^i'ipQE) h^s a somewhat asjTiimetric shape and 
a tail that extends towards positive values of ipQE- In contrast, the RFG (see Eq. (fT^ ) is 

^ In the figure and henceforth the QE scahng variable is denoted V'qb to distinguish it from the scaling 
variable used in the A region; see below. 
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symmetric when plotted as in the figure, is hmited strictly to the region —1 < iPqe < +1 
and has a maximum value of 3/4, while the empirical scaling function reaches only about 
0.6. 

One source for the difference could be typical mean-field dynamics in the initial and final 
nuclear states involved; however, for the kinematics of interest in the present work, both 



relativistic mean-field theory 



17 



20 



23| and relativized shell-model studies |l9j appear 



to provide only rather modest differences from the RFG predictions. As stated above, non- 
relativistic modeling is quite incapable for such kinematics (see, for instance, 

Another source for the differences seen between the RFG or mean-field descriptions and 
the empirically determined scaling function arises from high-momentum components in re- 
alistic wave functions which may be large enough to produce the results shown in the figure, 
although much more work, especially in a relativistic context, is required to put this on 
solid footing. For the present we limit our approach to phenomenology and take the scaling 
function from experiment. It should be emphasized at this point that much of the inability 
of typical modeling to account for the inclusive response at these kinematics appears to stem 
directly from the inability of that modeling to account for the results in the figure. We shall 
see below that when the empirical scaling function is used one obtains a good representa- 
tion of measured (e, e') cross sections and, therefore, that one's confidence in proceeding to 
predictions for neutrino reaction cross sections must be raised. 

Unseparated experimental quasielastic cross sections have been measured for several nu- 
clei over a large range of momentum transfer (for a compilation see 14]). These data have 
been used to determine the values oi kp for the nuclei considered (in p| a table with numer- 
ical values is given) and, since the evolution of kp with the nuclear mass number is slow, 
these values can easily be interpolated for any nucleus of interest. 

As was pointed out above, the description of the experimental scaling functions involves, 
besides the choice of the proper kp that sets the overall momentum scale, the use of an 
energy shift E shift that in an average way accounts for the nucleon removal energy. This 
small correction has been included in 0, 0] in the analysis of the data, and gives a table 
with the numerical values for various nuclei. 

We should add that, in order to obtain f'^^ from the cross sections, one had to divide 
by the single-nucleon cross sections obtained from e-p and e-n scattering. We have used the 
Hohler parameterization 8.2 j^^. In the range of q of interest here, this parameterization 
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agrees with more recent parameterizations fitted to a somewhat more extensive data set. 



B. Scaling in the region of the A peak 

Following the framework of 0, Is^, , let be the mass of a generic nucleon excitation 
and /i* = m^lm^\ hence /i* = 1 for quasielastic scattering and /i,, = ra/^lmj^ = /ia for 
electro-excitation in the A region. Introducing 



ft 



i ("2 - 1) 



= 1 + /5*/r 



(24) 
(25) 



we generalize the dimensionless scaling variable of the quasielastic peak as follows: 



^ (i^VpI + 1/^ - Ap* - 1 



V2 I +1 A > A" 
1 A < A? 



(26) 



which vanishes for 



A = A 



1 r 



or, in dimensionful variables, when 



(27) 



(28) 



When p* = 1 one recovers the QE answer in Eq. ()19p. where = 0, p* = 1 and at the 
peak lUq^ = ^/rnj^ + — m^r. As in the previous subsection where the QE scaling variable 
was discussed, here also we include the small energy shift E shift by making the replacement 
uj — E shift with A ^ A' and r ^ r' as before. Again these replacements are made 



UJ 



in the above equations to yield a generic shifted scaling variable ip'^, and specifically for use 
in the A region the shifted scaling variable ■j/'a- 

When considering the N ^ A transition structure functions we change notation from 
the general quantities p*, ip^^, etc., to /3a, Pa, "^a, etc., and in addition introduce 



7a = ^ (Pa 



and 



— \/r+{r + (3Ay 

A^A ^ 



(29) 



(30) 
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which allow us to define 



z/f = (1 + ^a) {t + 7a) 



(31) 



and 



r 



Then the iV — > A single-baryon responses will read: 



W2 l^j 



At 



/^A 



where the magnetic, electric and Coulomb form factors (following 

G'A/,A(r) = 2.97^7a(Q') 
Gea{^) = -0.03(7a(Q') 
GcAr) = -0.15G'MA(r), 



(32) 

(33) 
(34) 



are taken to be: 



with 



and 



GEp{r) 



1 



(35) 
(36) 
(37) 

(38) 
(39) 



[1 + 4.97r]^ ' 

namely the dipole parametrization of the proton (elastic) electric form factor (see also 
Sec. HvT^ . 

As for the quasielastic region, in the A domain we ignore terms of order rjp. In this 
approximation (as above, denoted by the subscript 0) the RFG longitudinal and transverse 
N A responses will read: 

R^{^, A)o = [(1 + rpl) w^{r) - (r)] x fnMi^^) (40) 



(41) 



where Aq is given in Eq. pSj). As usual, for electron scattering one should add the contri- 
bution obtained from Eqs. (j40|l and (j4T| computed with M = Z and the p A^ structure 
functions to the one where Eqs. (jinjl and pT|l are computed with J\f = N and the n — > A° 
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responses. Since these processes are purely isovector, clearly this is equivalent to using 
M = A with one set of the structure functions. 

Again, using the guidance provided by the RFG, the above procedure is easily generalized 
to the experimental response in the A region. Here, as long as density-dependent corrections 
{i.e., the corrections that go as rfp) are ignored as above, one should divide the experimental 
inclusive electro-excitation cross section in the A region by 

5^ = aM [vLGt + vtG^] , (42) 

to get the scaling function A). By comparing Eqs. ()4U|) and ()4H) with Eqs. ()13|) and 

(fTl]) one gets for the functions G^rp the expressions 

Gt = Yr W{{^ + ''pl)^2{r)-w^{r)]]+0{r^l) (43) 
G^ = -[^{wt{r)}]+0{vl). (44) 

Note that in Eq. (jlH) one has (r), whereas in Eq. (j221) the factor r in Wi{t) = tG\^{t) 
has been taken out in front. As before, one should take the results for reactions with protons 
with the appropriate form factors and with M = Z and add it to the results for reactions 
with neutrons again with the appropriate form factors, but now with M = N . Finally, to 
get the superscaling function /^{n, A) one multiplies F^{k, A) by kp. 

With the formalism in hand, we now proceed in a manner that is analogous to our treat- 
ment of the data in the QE region, however, now focusing on the A region. In order to 
isolate the contributions in the A region, we have subtracted from the total experimental 
cross sections (with Coulomb distortion effects incorporated) the quasielastic cross section 
recalculated using the universal f^^i^pQE) introduced above. That is, we remove the impul- 
sive longitudinal and transverse contributions that arise from elastic cN scattering, leaving 
(at least) MEC effects with their associated correlations and impulsive contributions aris- 
ing from inelastic eN scattering. As discussed earlier, the MEC effects will be ignored in 
the present work as they are believed to provide relatively small corrections, and thus this 
yields, at least for < 0, a response which is largely dominated by the A. For energy losses 
beyond the maximum of the A peak, other resonances and, at the larger values of q, the 
tail of deep inelastic scattering contribute. As a consequence of the different g-dependencies 
of the various contributions, it has not been possible using the present approach to analyze 
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this ip'^ > region further. The region of vahdity of /^("^Aa) therefore will be restricted to 

As shown by Eqs. (|33p and (|34j) . the determination of /'^(k, A) involves a division by a 
combination of the M-, E- and C-contributions with their appropriate g-dependence. For 
the latter, we employ the parameterizations given in Eqs. (|35H37j) used by Amaro et al. js^. 
Obviously, the main contribution is due to the Ml-term, the CO and E2 contributions to the 
cross sections being minor. 




-3-2-10 1 2 3 



FIG. 2: Averaged experimental values of /^(V'a) together with a phenomeno- 
logical fit (whose validity is restricted to < 0). 



In Fig. 121 we show the resulting f^(ip'^) extracted from the high-quality world data for 
inclusive electron scattering from ^^C and ^^O in the QE and A regions. These data span 
energies extending from 300 MeV to 4 GeV and scattering angles from 12 to 145 degrees, 
depending on the beam energy. For this determination the larger-angle and higher-g data 
are of particular importance. At small angles and lower q, the A contribution is small, and 
often not present in the available data due to limited coverage in energy-loss. As for f'^^, 
the experimental values of have been parameterized by a simple analytical function. We 
show in Fig.|21the averaged experimental values together with this fit. As pointed out above, 
the validity of the fit is restricted to cj- values below the peak, i.e., ip'^ < 0. 

The data appear to scale reasonably well up to the peak of the A, namely, the point 
where ■j/'a — 0' although clearly for still higher excitation energies the scaling is broken 
by processes that are not well represented via A-dominance. There is also some excess at 
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large negative ip'/^ which breaks the scahng to some degree — this is t 



lought to be due to 



contributions from MEC and their associated correlations jsiy, 10,|29, 30, 31, 32, 33, Is^, 
as discussed above. For reference we note that, for electrons of 1 GeV scattering from ^^C, 
the quasielastic peak (where — 0) occurs at = —1.8, —1.2 and —1.0 for 6 = 45, 
90 and 135 degrees, respectively. From Fig. we see that f^^ peaks at roughly 0.6 and 
thus these non-scaling contributions typically occur at the 10-15% level in the total cross 
section. Below we discuss our expectations for the uncertainties incurred for our predictions 
of charge-changing neutrino reactions when we ignore such effects (see Sec. fV]) . 

In passing it is important to note that in the present study we have simply taken the 
residual scaling function from experiment. While similar to f^^ it differs in detail: it 
is somewhat lower, is shifted slightly and is more spread out over a wider range of scaling 
variable. This is perhaps not unexpected, since implicit in this approach is the fact that the 
A brings with it its own width and shift. Only with a more microscopic model could one 
hope to be able to deconvolute these from the total response and see whether or not the 
underlying scaling function is indeed the basic f'^^ deduced above. Such an approach will 
be pursued in the future, although it only becomes practical when the MEC contributions 
are under control. For the present we limit the analysis to using two different functions f^^ 
and both deduced from phenomeno logical fits to electron scattering data. 

With the above ingredients, it is then possible to recalculate for every nucleus, incident 
electron energy and scattering angle the inclusive cross section for u below the maximum 
of the A contribution. In order to demonstrate this, in Figs. OHSlwe show the experimental 
responses together with the calculated response obtained using the parameterized /'^^ and 
In particular, we have studied the accuracy of the predicted response using (e, e') for 
^^C and ^^O and for a variety of momentum transfers, since these are the most relevant 
nuclei for the MiniBooNE and K2K/T2K neutrino oscillation measurements discussed in 
the introduction. For the data sets which do cover the A region, typical deviations are 10% 
or less. 
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Z = 6, A=12 

6=1500 MeV, 0=13.5° 
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cj (MeV) 

FIG. 3: Experimental (e, e') cross section for ^^C at an incident electron 
energy of 1.5 GeV and a scattering angle of 13.5 degrees, together with the 
calculated result obtained using /'^^ and f^. The dashed curve is the QE 
contribution and the solid curve is the total including the A. 



2.0 




200 400 600 800 

u (MeV) 

FIG. 4: As for Fig. |S1 except now for an electron energy of 1.3 GeV and a 
scattering angle of 37.5 degrees. 



IV. CHARGE-CHANGING NEUTRINO REACTION FORMALISM 

Among several options available (see for instance [20. 21. 

0,0, Q) 

we choose to write 

the charge- changing neutrino cross section in the target lab. frame in the form 



dndk' 



(45) 
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FIG. 5: Experimental (e, e') cross section as above, but now for ^^O at an 
incident electron energy of 0.88 GeV and a scattering angle of 32 degrees, 
together with the calculated result obtained using f'^^ and f^. 



where x = + for neutrino-induced reactions (for example, ui + n —>■ £ +p, where i = e, fi, r) 
and X = ~ for antineutrino-induced reactions (for example, Vl + p ^ + n). In Eq. ()45j) 

(46) 

where G = 1.16639x10"^ GeV~^ is the Fermi constant, 6c is the Cabibbo angle (cos^c = 
0.9741) and the generalized scattering angle 6 reads 

tan2^/2 = i^, (47) 

with 

v, = {e + e'f-q^ = Aee'-\Q^\ . (48) 

Henceforth we shall assume that m = = 0, but will always keep m' nonzero. 

One additional issue arises in computing the neutrino reaction cross sections having to 
do with the fact that the charged leptons in the final state are not plane waves, but are 
influenced by the Coulomb potential provided by the nucleus. This implies that the 4- 
momentum of the scattered lepton, K''^ = (e',k'), is the local quantity, in the sense that 
the 3-momentum k' and energy e' = V m'^ + k'^ are determined using the sequence of steps 
outlined in Sec. |nj culminating in Eqs. (fTUI) and pi|) for these variables. However, the 
asymptotic energy-momentum is not the same as the local quantity. Following standard 



[Lr COSfc 

2^ 



k' cos^/2 
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procedures (see, for instance, |2l|, |41|) the Coulomb interaction can be incorporated, at least 
approximately, by shifting from (e',k') to asymptotic energy-momentum {e'^,k'^): 



CO 

I 



D{k')k' 



where 



D{k') = l-x 



3Za 
2Rk' 



(49) 
(50) 

(51) 



and R = 1.2A^/^ is the effective charge radius of the nucleus being studied. Thus, our 
procedure is to calculate the cross sections using the kinematics as discussed in Sec. ITTland 
then at the end present the results in terms of the asymptotic energies and momenta obtained 
in this approximate manner. The only remaining issue is that the local calculations must 
also be multiplied by the density-of-states factor [D{k')]^^ to obtain the results we present 
in Sec. El 



The nuclear structure dependent quantity JF^ may be written as 



VccRcc + '^VclRcl + VllRll + VtRt 



X 



(52) 



that IS, clS cl generalized Rosenbluth decomposition having Charge-Charge (CC), Charge- 
Longitudinal (CL), Longitudinal-Longitudinal (LL) and two types of Transverse (T,T') re- 
sponses. Next we expand these response functions into their vector and axial-vector contri- 
butions according to 

(53) 
(54) 
(55) 
(56) 
(57) 



Rcc 


TjVV 

— ^cc 


+ Rcc 


RcL 


T3VV 

- ^CL 


+ R^ct 


Rll 


TjVV 

- ^LL 


+ Rit 


Rt 


dVV 

itrp 


+ R^^ 


Rfi 
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The lepton kinematical factors are the following: 



Vcc = I- tan^ e/2 ■ 6^ 

VcL = iy+-tanH/2-6^ 
P 

Vll = z/^ + tan^ 
1 



2u 



Vt 
Vt' 



p + tan^6'/2 



P 



) tan^ e/2 (^u + ^pp' -6^) -5^ 



P 



-ta.n^6/2 



(1 - vp' . 6') 



(58) 
(59) 

(60) 

(61) 

(62) 



Here rriiy has been assumed to be zero, and the entire lepton mass dependence occurs via 
the dimensionless parameter 



6 = 



m 



^/\Q' 



Moreover, in the above we have defined 



_ uj A 
q K 

m 

p ~ 



T 



P = 



e + e' 



(63) 

(64) 

(65) 
(66) 



which turn out to be related as follows: 



P 



tan 9/2 



p + tan'^e/2 



(67) 



and lie between zero and unity. In passing we note that, using the generalized scattering 
angle in Eq. (gTj), 



\Q^\ = 4ee'sin2^/2 
vo = 4ee'cos^^/2. 



(68) 
(69) 



Finally, for use later we define 



Vl = Vcc - 2iyVnL + ly'V, 



CL 



'LL- 



(70) 



For comparisons with finite-mass corrections to electron scattering see |4 
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In the extreme relativistic limit (ERL), namely m' 0, the kinematical factors are 
obtained by observing that in this situation 9 becomes 9 and that all of the terms containing 
5 can be dropped. One then gets 

Vcc Vcc = 1 (71) 

VcL ^ VCL = 1^ (72) 

Vll ^ VLL = (73) 

Vl ^ VL = (74) 

Vt ^ v^ = }-p + tan^ 9/2 (75) 



Vt> vt' = tan^/2y'p + tan^ 9/2 , (76) 

where the last three kinematical factors coincide with those employed in electron scattering 
(spin observables, coincidence electron scattering, parity- violating electron scattering, etc.). 

In the case of the /i = (C) and 3 (L) components of the vector current, which is 
assumed to be conserved, it is possible to collapse the contributions down to a single term: 
in particular, one has 

R^al = -yRVc (77) 
Rll = ^'Rll ■ (78) 

Since everything of the purely polar-vector type can be related to a single response, tra- 
ditionally we call this the longitudinal contribution, defined by the equation -R]^^ = B^xj- 
Expressing the sum of the (/i, i^) = (0,0), (0,3), (3,0) and (3,3) contributions, one then 
ends up with a single term 

VccRi:'c + ^VclRH + V,,Rll = V,RY ^ ■ (79) 

This collapse into a single expression does not occur for the AA terms. Indeed there one 
has 

yccRcc + 2V"cL-RcL + ^llRll = ^cfh ■ (80) 

To complete the analysis one should add the two contributions (//, i^) = (1,1) and (2,2), 
which yield 

% + ^T^] = Xt , (81) 
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and, as well, consider the V/A interference term where (/i, i^) = (1, 2) and (2, 1) 
The full response will then be (see Eq. 



(82) 



(83) 



A. Single-nucleon responses in the QE region 



The single-nucleon responses in the QE region all begin with the basic vector and axial- 
vector currents involving N ^ N matrix elements: 



fv 



3a 



u{P') 
u{P') 



2m N 
1 



2m 



-GpQ^ 



N 



U{P), 



(84) 
(85) 



with = P'^ — P^. Indeed, in the scaling analyses discussed in Sec. 1111 A| the former 
was used together with the usual relationship between the Dirac/Pauli form factors and the 
Sachs form factors, Ge = Fi — TF2 and Gm = Pi + to obtain expressions such as those 



in Eqs. fjl6ll7p . The total current is then = jy — j^. 



In fact, for the purely polar- vector contributions we have 

1 r ,-..12 
P L 

2r 



r>VV 

JrCrp 



(86) 
(87) 



which were used above. Here G^^\,j are the nucleon's EM isovector form factors. We re- 
emphasize that, as written, these results contain effects from the motion of the nucleons 
in the nucleus up to first-order in rip and only terms of order rjp and beyond have been 
neglected. 

For the axial-vector current proceeding from Eq. ()85p one has 

2 



R, 
R 



cc 

AA 
CL 

AA 
LL 



P 



G 



'(1) 



P 



G 



./(I) 



G 



'(1) 



oAA 
rCrp 



1 

P L 

2fl + r 



G 



(1) 



(88) 
(89) 
(90) 
(91) 
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defining the following combination of axial-vector and induced pseudoscalar form factors 

GT{r) = G^^{T)-rG^^{r), (92) 

namely, the axial-vector analog of the relationship between the Dirac, Pauli and Sachs form 
factors. It is of importance to realize that in the ERL (see Eqs. (j71H73|) ) X^^^ — > 0: hence 
this term crucially depends upon the final lepton mass, yielding for the C/L single- nucleon 
QE response 

■ (1 + 5^)52, (93) 
Finally the V/A interference term is 

E^f^ = 2v/^(rT7)G£)G(J) . (94) 



X^f^ =tan^^/2 



G 



7(1) 



B. Single-nucleon responses in the A region 

In this subsection we consider the elementary reactions 



















"AO 






-A" 



The associated currents are 
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(95) 
(96) 
(97) 
(98) 

(99) 



where the isospin factor T is for A^'^ and A" production and 1 for A+ and A° pro- 
duction, u[f^\p',s') and u{p,s) are the Rarita-Schwinger and Dirac spinors for a A and a 
nucleon with momenta p' = p + q and p and spin s' and s, respectively. For the vertex tensor 
we take 12211 



- r'V 



+ 



'cl 

TUN 



(9^ 



- q"!^) + -4- {9'^^q ■ v' - 9 V^) + ^ {9"'^(i ■ V - q> 



75 



m 



N 



m 



N 



(100) 



25 



We recall that CVC implies = and PCAC yields = C^{fil + 4t)-\ with /i^ = 
mT^/niN, m^r being the pion mass. 
The hadronic tensor 

w'^'' = fLATT{J'^^{q)r{q)}, (101) 
with /iA = Tn/^/rriN as above, can be rewritten in the form 

w^'' = T'^^aTt {P^^ip') (7ort"'^7o) A{p)T^''} , (102) 

where 

s' 

= ^ 9f^a - + q7/37a (104) 

is the Rarita-Schwinger projector, while as usual 



2mAr 

is the nucleon projector. 

A lengthy calculation then yields 



A(P) = ^-P^ (105) 



where 



w^'' = w^''y + w''^^ + w^\ (106) 



<v = -wivig^'' + ) + W2v{v^ + Pi^^^niv"" + P^^^n , (107) 

r 

T 

-UiA^ + U2A(ft:V + ^'''«") (108) 

T 



and 



= 2^w;3e"'^^^r/„«:^ (109) 

are the Vector- Vector, Axial-Axial and Vector-Axial interference hadronic tensors, respec- 
tively. The functions Uia and U2a reflect the non-conservation of the axial-vector cur- 
rent and will be discussed in the Appendix. The functions Wi entering above are ob- 
tained by performing the traces in Eq. ()101|1 and using the on-shell conditions 77^ = 1 
and T) • K = rpA = t + ^(/i^ ~ 1)- The results are collected in the Appendix. 
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C. Nuclear cross sections and response functions 



To compute the cross section in Eq. the factor J-"^ in Eq. ()83|) is required, and thus 
the nuclear response functions Ri are needed. The last can be expressed in terms of the 
nuclear tensor W'^'^ according to 

RV = (110) 

R^^ = WZ (111) 

Ret = -WTa (112) 

Rtt = WZ (113) 

RV = Kv + (114) 

Rt'' = W\\ + WI\ (115) 

R^;^ = -iW^\ . (116) 

Since in the QE domain the nuclear tensor is well-known, here we focus on the A sector, 
where using Eq. (fT5|) W^^ reads 3^ 



= ^AoW(v^A)t/^^ (117) 

where M = Z for reactions on protons, (j95p and (j97|) . whereas M = N for reactions on 
neutrons, (|96|) and (|98|) . Again, only isovector form factors enter. 

As discussed previously (see Sec. IIIIBj) . the REG superscaling function in the A domain 
is given by 

W^A) = ^(i-v^i)^(i-^i) (118) 



and 

1 

where 



V^"" = — - — \^ w^\e)de, (119) 
ep — eo 



eo = ^^\jp\ + 1/r - ApA (120) 

represents the minimum energy a nucleon should have in order to take part into the A 
electro-excitation process and ep = ^yl + rjp is the Fermi energy. Note that, from Eq. 
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one has 



7- (eo - 1) 



1/2 



X 



+1 \>\\ 

-1 A < 



(121) 



with \\ = \ 



^JlJi\ + Ah? - 1 using Eq. (jTFj) . 
The required components of the U^^ tensor can be exphcitly computed using Eqs. ()107t - 
and performing the integral in Eq. ()117j) . They turn out to be 

— [-wiv{t) + (1 + tp\)w2v{t) + V{k, t)w2v{t)] 

T 

(122) 



rrOO 



rrOO 



T 

A^ 



[-u;iA(r) + (1 + rp^^)w2A{r) + P(k, r)ii;2A(r)] 



UiAir) + A(eF + eo)u2A{r) 



(123) 



rr03 



-] { — [-WiAir) + (1 + rp^)w2A(r) + V{k, r)w2A(r 
A2 + ft:2 



Mia(t 

r 



2A 



-(eiT' + eo) - TpA 



U2A{r) 



(124) 



7-r33 



A\' 



[-w^ia(^) + (1 + rp\)w2A{r) + ^^(fi;, r)tf;2A(^ 



rA2 



UiA{r) 



A2 



[A(ei;' + eo) - 2tpa] «2A(r) 



(125) 



for pieces of the tensor having p oi v equal to or 3, while for transverse projections one 
has 



Uyy + Uyy 



^AA + ^1 A 

TTl2 

•^VA 



2wiy{T) + P(k, t)w2v{t) 

2wiA{r) + V{k, T)w2A{r) 



2zJt{1 + tpI)[1 + V'{k,t)]w3{t) 



(126) 
(127) 
(128) 



where 



V{h,t) 



T 
t? 
1 
K 



- (e| + eoCF + el) + ApA (cf + eo) + >? p\ 



ApA + - (cf + eo) 



1 . 



(129) 
(130) 
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Thus we are now in a position to assemble the various factors and provide predictions 
for the charge-changing neutrino reaction cross sections. The single-nucleon cross sections 
in both the QE and A regions are given above and, from the scahng analysis presented 
in Sec. 11111 we have experimentally determined scaling functions f'^^i'ipQE) /^(V^a)- 
As usual, we have used scaling variables "ipQE and in which we have included a small 
energy shift by replacing u with u — Eshift = u', \ ^ A', etc. For charge- changing neutrino 
reactions one proceeds to analogous states in the neighboring nuclei, which are shifted from 
their positions in the target nucleus at least by the Coulomb energy. Accordingly, while a 
very small effect, in this work we have modified Eshift from its value as obtained in studies 
of electron scattering by these additional amounts. In the case of mass-12 this means that 
upon comparing the excitation energy of the analog of the nitrogen and boron ground states 
{i.e., the energy 15.110 MeV of the lowest J'^T = state in carbon), we should add 
16.827 — 15.110 = 1.717 MeV when making transitions to nitrogen (neutrino reactions) 
and subtract 15.110 — 13.880 = 1.230 MeV when making transitions to boron (antineutrino 
reactions) from the canonical value of EsUft = 20 MeV. These are then used in the definitions 
of ip' in each case. Note that we use the correct ground state masses of the three nuclei in 
establishing the lepton scattering kinematics discussed in Sec. m 

Finally, in computing the results to be presented in the next section the following form 

'actors have been employed. For the vector sector we use the Hohler parameterization 8.2 

3fi| and for the axial-vector sector we use 

GUr) = (1 + A^r)-^ (131) 
G^}\t) = saGUt) (132) 
G^hr) = Tn7—zG^}\r) (133) 

with A;^ = 3.32 (corresponding to axial-vector mass Ma = 1032 MeV), A^ = {2mN /mT,Y = 
180 and gA = 1.26. 



V. NEUTRINO CROSS SECTION PREDICTIONS 

The inclusive charge-changing neutrino reaction cross sections ^^0(1/^, /i~) and ^^C(z/^, /i"*") 
that result from the scaling analyses presented above are shown in Figs. Inland [71 respectively. 
A neutrino or antineutrino energy of 1 GeV has been selected as being representative of 
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FIG. 6: Charge-changing neutrino reactions on C for 1 GeV neu- 

trinos and neutrino-muon scattering angles of 45, 90, 135 and 180 degrees. 
The cross sections are plotted versus the final-state muon momentum k' . The 
dash-dotted curves give the QE contribution, the dashed curves the A con- 
tribution and the solid curves the total. As discussed in the text, results for 
values of k' lying below the A peak (higher excitation energies than that of 
the A) must be viewed with caution. 



kinematics where the scaling approach should be expected to work well and a selection of 
scattering angles (between incident neutrino or antineutrino and produced charged muon) 
has been made. Note that, since the predictions shown here are given as functions of 
the muon momentum /c', the QE peak lies to the right {i.e., lower excitation energy) of 
the A peak (higher excitation energy). Note that the Coulomb distortion correction from 
Eq. has been included here, although for simplicity we have labelled the figures using 
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FIG. 7: As for Fig. |H1 but now for antineutrino reactions {v^, fi~^ 



k' rather than k'^. As discussed above, the predictions at momenta to the left of the A 
peak (excitations lying above the A region) are unreliable, since our scaling approach does 
not fully account for meson production, including resonances other than the A, and DIS 
processes. 

Corresponding angular distributions are shown in Fig. |H1 for kinematics chosen to lie at 
the peaks of the QE (solid curves) and A (dashed curves). 

One striking feature of the results is the dramatic differences seen in comparing neutrino 
and antineutrino cross sections at backward angles. The latter are typically about two orders 
of magnitude smaller than the former under those conditions. This is due to the fact that 
the transverse (vector and axial-vector) contribution to both the QE and A responses is 
accidentally roughly the same in magnitude as the one arising from the V/A interference for 
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FIG. 8: Angular distributions for the results in Figs. El and [7| at the tops of 



the QE and A peaks. Since the neutrino or antineutrino energy is fixed to 1 



GeV and the kinematics are chosen to be at ip'r 



QE - 



and V-'a ~ 0' muon 



momentum k' 



must vary with 0. 



large scattering angles at the chosen energy. Specifically, for instance in the QE case, one 
has Xt roughly equal to Xt' (see Eqs. dHU), (IH21), (IHZ|), ^ and For neutrinos these 

constructively interfere, whereas for antineutrinos they tend to cancel and produce much 
reduced cross sections. Indeed, the cancellation is so severe that the VV(CC, CL and LL) 
and AA(CC, CL and LL) terms can yield significant contributions to the total cross section. 
In the QE region these VV terms yield as much as 1/3 of the cross section, while the AA 
terms are negligible. In contrast, in the A region the reverse is true, with the AA(CC, CL 
and LL) terms even providing the majority of the cross section at large angles in this case. 
Small changes in the model (for example, the inclusion of MEC) could have very large effects 
on the predictions for antineutrinos and hence the results shown in this case, especially for 
large scattering angles, should be viewed with caution — the important observation is that 
the antineutrino cross sections are predicted to be strongly suppressed due to the accidental 
cancellation. 

For neutrinos, where no such strong cancellation occurs, the cross sections are typically 
dominated by the VV(T), AA(T) and VA(T') contributions. In the QE region the VV(CC, 
CL and LL) pieces contribute only about 5% at 45 degrees and fall to negligible corrections at 
backward angles, whereas the AA(CC, CL and LL) contributions are negligible for all angles 
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considered. In the A regions the converse is true: the AA(CC, CL and LL) pieces contribute 
about 5% or less at 45 degrees and fall to negligible corrections at backward angles, whereas 
the VV(CC, CL and LL) contributions are negligible for all angles considered. 

The fact that the underlying transverse vector and axial- vector matrix elements are com- 
parable in magnitude for the kinematics being explored in this work (hence the antineutrino 
suppression discussed above) has consequences for the uncertainties expected in the pre- 
dictions made here for neutrinos. We recall from our treatment above that effects from 
MEC and their associated correlations are not present for these predictions, since they were 
ignored in analyzing the electron scattering cross sections. Note, however, that the error 
incurred by this may be less here than for electron scattering: the transverse contributions 
to the neutrino reaction cross sections involve both polar- and axial-vector matrix elements 
with both one-body, impulsive contributions (included here) and two-body MEC/correlation 
contributions (not included). From past studies one knows that the latter occur primarily 
in the transverse channel, but not in the longitudinal channel for electron scattering at the 
kinematics of interest here. That is, the vector MEC/correlation effects occur primarily in 
the transverse channel for QE and A kinematics at high energies. In contrast, due to the 
factor 75 that enters in axial-vector currents, the converse is true for axial-vector contribu- 
tions — accordingly, to leading order one does not expect large corrections of this type for 
the axial- vector contributions. As a consequence, the residual effect seen in Fig. |21 at large 
negative ip'^ and attributed to MEC and their associated correlations which amounted to 
roughly 10-15% of the total cross section in that region of inelasticity measures the uncer- 
tainty in the vector contributions. Since these are only roughly 1/2 of the total for neutrinos, 
with the axial-vector transverse contributions accounting for the other half, the overall im- 
pact of the neglected MEC/correlation contributions is likewise only half as large, namely, 
providing less than 10% uncertainty to the neutrino cross section predictions made here. 

For completeness in Fig. 121 we show a comparison of a typical neutrino reaction cross 
section obtained using the full results with f^^ and deduced directly from electron scat- 
tering data with the cross section obtained using the RFC. As noted earlier, modeling done 
on the basis of mean- field theory is very similar to the RFC result, since at the relatively 
high energies of interest the dynamical effects embodied in an effective mass are expected 
to be small (most recent relativistic mean-field theory studies predict that m* reverts al- 
most to the nucleon mass for the kinematics of interest). Likewise, relativized shell model 
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FIG. 9: Neutrino reaction cross sections as in Fig.^lfor ^ = 45 degrees, showing 
a comparison of the full results obtained using the empirical scaling functions 
f'^^ and obtained as discussed in the text with results obtained using the 
RFG scaling function Jrfg (heavier lines). The former lie somewhat lower 
and extend over a wider range in k' than the latter. 

predictions are close to the RFG predictions, and, moreover, RPA correlations are expected 
to be relatively small for the high energies involved. Thus, the RFG predictions effectively 
represent a larger set of models. As can be seen in the figure, all therefore differ significantly 
from the scaling predictions. Given the success of the scaling approach in studies of inclusive 
electron scattering for the kinematic region under study we expect that neutrino reaction 
cross sections also obtained using scaling ideas to be more robust than those based directly 
on existing models. 

VI. CONCLUSIONS 

We begin by summarizing the approach followed in the present study. 

The first step has been to explore the scaling behavior of inclusive electron scattering 
for relatively high energies (several hundred MeV to a few GeV) in the kinematic region 
extending from the scaling region which lies below the QE peak, through the QE peak 
region and up to the peak where A-excitation is the dominant process. Upon examining the 
longitudinal contribution one finds superscaling, i.e., it is possible to find a scaling function 
jQE ^Jiich, when plotted versus an appropriate scaling variable ipQ^, is seen to superscale. 
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This means that the results are found to be independent both of the momentum transfer 
(scahng of the first kind) and of the particular nuclear species (scaling of the second kind). 
The assumption is made that this universal scaling function embodies the basic nuclear 
dynamics in the problem. Implicit in this, and borne out by modeling, is the observation 
that apparently MEC effects with their associated correlations and inelastic processes such 
as A excitation are not large contributions for the longitudinal response. 

The second step has been to use the universal scaling function f^^ to obtain that part 
of the transverse response that is due to impulsive processes, namely, those which arise only 
from elastic eA^ scattering from nucleons in the nucleus. Upon subtracting this contribution 
from the total transverse response one finds a residual which is assumed to be due to the 
MEC/correlation effects and to inelastic eA^ scattering from nucleons in the nucleus. From 
modeling we know that such effects are predominantly transverse and hence naturally occur 
in this channel, but are not very important in the longitudinal channel. Moreover, again 
from modeling of the various processes, we expect that the MEC/correlation effects are 
relatively small corrections for the kinematics of interest and accordingly attribute most of 
the residual to impulsive, inelastic eN scattering, especially to contributions which arise 
from ^ A transitions. When the residual is analyzed in terms of an appropriate scaling 
function by dividing the residual response by the elementary — * A cross section and 
plotted versus an appropriate scaling variable which incorporates the kinematics of the 
inelastic transition we again find a reasonably successful new kind of scaling, at least for 
kinematics where the concept is expected to work. Specifically, this success is only found 
for excitations up to the peak of the region where the A dominates, but not beyond: this 
is not unexpected, since the approach taken in the present work has been tailored only to 
work when the A provides the basic driving process. Also, deviations are seen in the region 
where the QE and A responses overlap and there one does expect corrections from MEC and 
their associated correlations, while small, nevertheless to be necessary for a fully successful 
representation of the total response. A check is made by reassembling the complete inclusive 
(e, e') cross section for the kinematic region of interest, typically finding errors of 10% or 
less. 

Thus, the first goal of this work has been met: we have found a very good representation of 
inclusive electron scattering at relatively high energies for the region of excitation extending 
up to the peak of the A. It should be noted that direct modeling [i.e., relativistic modeling, 
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since non-relativistic approaches are known to fail badly for the kinematics of interest) yields 
electron scattering cross sections that are valid only at about the 25% level or worse, while 
here our goal has been to use the scaling approach to do better. 

The second major objective in this study has been to predict charge-changing neutrino 
and antineutrino cross sections for the same range of kinematics. Using the scaling approach 
in reverse, we take the empirically determined scaling functions f^^ and together with the 
appropriate N ^ N and N ^ A charge- changing weak interaction cross sections to obtain 
the inclusive uA and uA cross sections for the case of ^^C. Given the above statements 
of where the approach taken in the present work should be valid we believe that these 
predictions should be the best currently available for few GeV neutrino reactions in the 
kinematic region that includes the full QE response and the A response up to its maximum. 
From our analyses of MEG contributions and how they enter in the relevant vector and axial- 
vector responses, we expect that corrections from such processes account for only about 10% 
of the total cross section. Note that, while our focus has been on the case of carbon, we 
know from our previous studies of scaling of the second kind that it is straightforward to 
produce predictions for other nuclei as well. In particular, while only selected predictions 
are given in the present work, further results at different kinematics and for other nuclei 
may be obtained by contacting our collaboration [4^ . 

Finally, we mention our intentions for future work. Most straightforward will be to extend 
the scaling approach to obtain predictions for neutral current neutrino and antineutrino 
scattering from nuclei for similar kinematics; studies of this type are currently in progress. 
Also, given recent work on inelastic eA^ processes and their role in eA inclusive cross sections 
undertaken by the same collaboration, our intent is to explore a more microscopic model for 
all of these reactions. While this project is relatively straightforward as well, what is lacking 
before it can be realized is a completion of our on-going studies of MEG effects, together 
with the correlations they require due to gauge invariance; studies of this type are also in 
progress. Finally, there is the issue of explaining the specific nature of the scaling functions 
themselves. While there are indications that contributions from high Fourier components in 
the nuclear ground state are probably responsible for the detailed nature of these functions, 
a fully satisfactory relativistic treatment of them is presently lacking. Until one becomes 
available it appears that the best approach is to take the scaling functions directly from 
experiment, as we have done in the present work. 
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Appendix 



Performing the traces in Eq. (llOlj) one gets: 



Wiv 



^ ) ITT [4r + (/iA-l)' 



3 ; 



16 + 8r(/iA + 1) + (/iA + 1) (3/ii + 1) 



CP 



(l + AT-fil)Cl + il-AT-fj.l)Cl 



/iAC3^(l+4r-2/iA-3/ii) 



il + AT-ixl)C^ + il-AT-fxl)Cl 



(134) 



= {^]^{il+^r + 3fil)Cr 



+ 



4r 

3 / /iA 

4r + (/iA-l)' 



(1 + 4 r - 4/iA + 3/ii)(Cr + Cs^) + 8 r C, 



for the vector contributions, whereas they yield 

'r2\ 1 



3 ; 4/ii 



[4r+(/iA + l)1' 



(135) 



16 - 8r(/iA - 1) + (/^a - l)l3/xi + 1) 



Cf 



+ /^A 



il + AT-fil)Ct-2C^ 



+ /xaC3^(1 + 4r + 2/iA - 3/i^) 



(l+4r-/ii)Cf -2^5^ 



(136) 



T2\ 4r 

3 ; /ii 



w^2A = \ — ]—<,{l + ^r + 3fil)Cf 



+ 



4r+(/iA + l)' 



+ /iA 



(1 + 4r + 4/iA + 3/ii)C4"^ + 2C, 



(137) 



38 



UiA 



[4r + (/iA + l)^](C5^-4rC, 



3 ) IQt ii\ 

48 r2-(/i^-l)2 + 8r (/i^ + l) 



+ 4 r 



4/iA {Ci + ^l^Ct) (4r + /ii-l) 



- C6^[16r2 + (/ii-lf + 8r(/ii + l)] 



(138) 



^2- 

3"y 4r/ii 



«2A = 1 — )^^[4r + (/iA + in(C5^-4rC6^ 



X [8 r /iA (Cg^ + /iA ) + Ct (4 r - /ii + 1)] 



(139) 



for the axial-vector contributions and 



W3 



2C5^-(4r-/ii + i) a 



X 



(1 + 4 r + 4/iA + 3/ii) CI - /iA[(4r - /il + 1) - (4r + /il - 1) 



2(1 + 8r + 16r2 + - 3/ii) 



+ /iA(l + 4r - 4/iA + 3/ii)[(4r - /il + 1) Q - (4r + /il - 1) C, 



(140) 



for the interference pieces. Since the axial-vector current is not conserved, Wi and W2 are 
not sufficient to set up the AA hadronic tensor: hence two extra functions Uia and U2A1 
which vanish if = 4:tCq, come out from the traces. 
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For the empirical functions entering above we take 22 



Ctir) 



Ct{r) 



(1 + |Q2|/o.54GeV^2)2 






-0.3 
1.2 ( 1 



1.21IQ2 



i.2i|g^| 

2GeV^ + 



(1.28)2GeV2 y 

.21 \ -2 



1 + 



(1.28)2Gey2 



^^^■^m2 + |Q2| 4^ + ^2 
which, inserted into Eqs. p35lll40p . lead to 



Wiv 



W2V 



(^)J[4.+(..+im 



2i,^y2 



for the vector contributions 
'T2\ 1 



W^2yl 



3/4 



[4r + (^A + l)2] [(4r-^i + l)C^-2Q^ 



T2\ 1 



[4r + (/iA + in(C5^-4rC, 



X 



3 y 16r/i^ 
48r2 - (/xi - If + 8r(/.i + 1) 



+ 4r 



4/ii (4r + l)Ct - [16r2 + (^^ - 1)2 + 8r (/i^^ + l)]^^^ 



M2A 



(141) 

(142) 
(143) 
(144) 

(145) 
(146) 
(147) 



(148) 
(149) 



(150) 
(151) 



(152) 



4r/i2 



[4r + (/.A + If] (Ct - 4r C^) [8r^iQ^ + (4 r - /^i + 1)^ 



(153) 



for the axial-vector contributions, and 



W3 



T2\ Cl_ 



[4r + (/iA + 1)'] [2C^ - (4r - /i^^ + l)Ci 



(154) 
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for the V/A interference. 

One finds that uia and U2a (which arise from PCAC) are neghgible, whereas the other 
functions are all significant. The latter are seen to fall strongly with increasing q. 
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